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ABSTRACT 


The eleven actuation system for a modern fighter 
airplane is modelled and analysed. Effect of various system 
parameters on the d’/namic response of the electro hydraulic 
actuator is studied. Impedance characteristics of the system 
are determined and their relation to the overall transfer 
function of the eievon actuation systesn is obtained. The model 
developed is integrated into the longitudinal flight control 
loop to investigate its effect on the airplane short period 
dynami cs . 

The response of the E H actuator under load conditions 
is oscillatory and lowly damped. This is far from the ideal 
first-order -lag response. However, for an optimised actuator 
position feedback the airplane normal acceleration response 
is good. though the pitch rate response suffers in 


comparison. 
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Chapter 1. INTRODUCTIOM 


1 • 1 I HTRODUCTI ON 

During the past decade, as the control pover 
requirements increased and the control reaction time 
decreased, more and more reliance was placed on fully powered 
hydraulic actuators. Electro hydraulic servo mechanisms are 
used to provide power boost, stability augmentation control 
inputs, or primary control through main and auxiliary control 
surfaces of the aircraft. These servos are capable of 
performance superior to that of any other type of servo. 
Large inertia and torque loads can be handled with high 
accuracy and very rapid response [13. 

A modern fighter airplane having a delta wing planform 
usually has a single set of control surface used for both 
pitch as well as roll control. These are known as .elevens 
C elevator aileron!!), in a typical airplane of this type* 
each eleven on either wing is split into two. Thus there are 
four elevens ~ two inboard and two outboard elevens. Each 
elevon is controlled by four servo actuators. In a fly by- 
wire aircraft the conventional mechanical link between the 
pilot stick and the actuators is eliminated and replaced by a 
multiplex flight control computer. Pilot commands are 
translated into electrical voltages which are fed to the E H 


ser VOS . 



The present work deals with the modelling and analysis 
of an elevon actuation system. A typical working electro 
hydraulic actuator on a modern fighter aircraft is considered 
for illustration. The actuator is of a linear dual tandem 
type controlled by a direct drive valve. 

1*2 LITERATURE SURVEY 

Literature on electro hydraulic systems is extensive. 
Text books by Merritt Cl 3 and Viersma [23 are excellent 
reference materials. Both these books present a rational and 
well balanced treatment of hydraulic components and systems. 
Hydraulic control is analysed by conventional linear analysis 
methods and general criteria applicable to “good** design of E 
H servos is evolved. In a nutshell* hydraulic servo design 


proceeds 

by first selecting 

the 

power element to 

meet 

the 

prime considerations of 

load 

and 

response 

capabilities. 

The 

r emai nder 

of the 1 oop 

is then 

sized to 

obtai n 

the best 

possible 

servo performance 

Ci . e 

. , accuracy and 

speed 

of 


response!) under the li nutations imposed by the power element 
dynamics. Noise rejection and stability alwa>rs limit the 
system bandwidth- In fact* it is desirable to keep the 
bandwidth at a nunimum consistent with the specifications. A 
reduced bandwidth usually simplifies compensation and, as 
peak power outputs are associated with high frequencies, it 
r el axes r equi r ement s on i ndi vl dual el ements 1 1 3 . 

McRuer £33 has analysed the actuator response to a 
spectral input to determine the general requirements for such 
a system. He concludes that ‘the actuator bandwidth should be 
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as low as possible^ consistent with suppression or actuator 
loop disturbances * effective linear isation of non linearities 
and outer loop stability and dynainic performance* of the 
aircraft. 

Flight control specialists are nowada^rs accepting high 
reliability in lieu of redundancy/ to increase the mean time 
between failures and decrease the complexity of the flight 
control system . Direct drive valves CDD\0 are a natural 
choise. Their use is based on the preniise that their 
reliability is so high that redundancy is unnecessaryC 43 , 
DDVs are being increasingly used instead of Electro Hydraulic 
Valves CEH\0 in ail modern fighter aircrafts. Hot much 
literature Is however available. Apart from some SAE papers 
C4], product information brochures of control system 
component manuf actur ers £53 are the only technical source 
available. Reference £43 and £53 briefly explain the working 
of the DDV and its advantages over the EHV. Reference £43 
also discusses other advanced concepts* like variable 
displacement hydi*"aulic control, electro mechanical servo 
systems, '•smart** valves etc. , which are yet in the 
experimental stage but nevertheless hold rich promise. 

Edwards £63 in a NASA technical note has analysed an 
electro hydraulic aircraft rudder servo. The transfer 
functions are developed and their approximate literal factors 
are presented. Aerodynamic and horn stiffnesses are taken 
into account. However the actuator is assumed to be rigidly 
mounted in the airframe. A noniineari ty, such as Coulomb 
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friction, presents difficulty in linear analysis of a system. 
EdwardsC6] has developed an equivalent linear model s^fherein 
an equivalent viscous friction coefficient has been 
considered instead of Coulon±) friction and compared it with 
non linear test results. However, no criteria have been 
proffered to determine such an equivalent viscous friction 
coefficient. 

Aeroser voelastic studies are now carried out during the 
design stage of all modern fighter aircraft to explore the 
possibility of active flutter suppression. A "real" actuation 
system is a necessary pai't of an Aeroser voelastic model. An 
actuator with finite impedance characteristics ix'istead of an 
infinite impedance actuator is therefore considered for such 
analysis C73. Benun CS3 has analysed a simplistic actuator 
model Cno masses; no dampers!) to illustrate the frequency 
dependant characteristics of the hydraulic actuator. 

1.3 PREVIEW OF THE WORK 

Chapter 2 describes the various elements of the electro 
hydraulic actuator — eleven system. The equations of motion, 
the trai'isfer function and the block diagrams are presented. 
Effect of change in values of the parameters of E H actuator 
on the system response is investigated. The instability of 
the load resonance poles due to eleven inertia is brought out 
and the dynamic pressure feedback is designed. 

Chapter 3 presents the importance of the impedance 
characteristics of the actuator - eleven system. A ph 2 /sical 
explanation is attempted for the actuator hydraulic 
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stiffiiess. The system impedance is derived and its relation 
with the transfer function between elevon deflection and 
command inputs of the elevon actuation system brought out. 

Chapter 4 describes the effect of putting the servo 
system in the overall flight control loop. The aim is not to 
design a control law, rather to see how the realistic model 
of the actuator developed interacts with the control law. 

Chapter 5 presents the conclusions drawn from the study 
and suggests possible extension of the present work.. 


CHAPTER a ELEVOM ACTUATION SYSTEM 


For t.he purpose of iirtalysis the dynamics of the 
actuator- -el evon servo system may be subdivided into: 
i:5 The direct drive valve dynamics 
ii’j The hydraulic actuator dynamics 
± i i 5 The actuator 1 oad dynami es 
A block diagram of a typical actuator - elevon servo system 
used in the flight control system of a modern fighter 
AirerAft ia ahown in Figur^r 2. 1 C13.* 

2a 1 DIRECT DRIVE VALVE 
2.1.1 I MTRDDUCTI OW 

Convent,! ona.1 2 — aLage el ecLro— hydraulic servovalves 

CEHVI> Are now being increasingly replaced by Direct, E)rive 
Val v^^sC DDV 3. EHVs achieve high power gain Lhrough hydraulic 
amplif'icaLion of* a small electrically generated torque. The 
operation oT the torque motor and hydraulic amplif’ier 
involves small Torces and displacements that makes the system 
vulnerable to jamming by possible contamination of" hydraulic 
fluid. Redundancy by multiplexing is therefore introduced 
which results in somewhat bulky and mechanically complex 
servo actuator systems. In addition» hydraulic system 
prc-saures are limited to around 4-SOOO psi as the EHVs are 
subject to leakage beyond this pressure . Development of the 
DDV seeks to overcome these limitations by means of increased 


mec hani cal r el 1 abi 1 i t y and eonser vat i ve redundant coi 1 
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conf igaration [4,53. The DDV has Taciiitated conTigurations 
having simplified control surface actuation systems and 
specification of higher hydraulic pressures C 8000 psi D as 
the construction of DDV tolerates pressures of this magnitude 
without leakage. 

The DDV, which is shown in Figure 8.2 , consists of a 

small high output electromagnetic torque motor which is 
directly coupled to a hydraulic metering arrangement C the 
main control vaiveD. A variable gap torque motor consists of 
two RECO C Rare Earth CObalt D magnets which set up a 
polarizing flux in the two air gaps. At null the flux in each 
air gap is balanced and there is no output. Control is 
effected by applying a comn\and signal to the coilCsZ) Casually 
there is one coil per multiplex channels wound around the 
outside of the motor. This current creates a control flux 
imbalance between the air gap and consequently a net torque 
on the motor armature. The motor is spring centered and 
therefore the opening of the valve is proportional to the 
command signal C4,5]. 

TRANSFER FUNCTION OF THE DDV 

The block diagram and associated transfer functions of 
a typical DDV used in Flight Control systems is given in 
Figure 8. 3. The rotary torque from the motor is coupled to a 
rotary control valve. For optimum performance a position loop 
is closed around the valve. This is also used to detect a Jam 
if it occurs and to generate additional voltages to chip away 
the jamming particle [43. This is a distinguishing feature of 
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the DDV. The 5^^ order transfer function which results from 
the block diagram is given in Appendix A CA1>. It may be seen 
from the values of the coefficients of higher order terms of 
s in the denominator of the transfer function that the 
transfer function may be approximated to a first order lag 
thus: 

ecs:) Ki 

= ca.iD 

V CsD Ts + 1 

o 

where, Kl = 0.0254 rad/voit 
T = 0. 0042 s 

A plot of the frequency response of the 5^^ order system and 

s t.^ 

the 1 order approximation of the DDV is shown in Figure 2.4 

s t 

The figure shows a good match between the 1 order and the 

5^^ order system up to 400 radxs. The range of characteristic 

frequencies of the eleven servo system is usually up to 10 Hz 

s t 

C 65 rad/sD. The 1 order system is therefore used to 
represent the DDV for further anaiysis.lt may be noted from 
the frequency response curves in Figure 2. 4 that the 
bandwidth of the DDV is 240 rad/s. 

2.2 HYDRAULIC ACTUATOR 
2-2.1 ACTUATOR DYNAMICS 

The hydraulic actuator considered is a dual — tandem — 
actuator. Figure 2.5 shows the schematic diagram of a typical 
valve actuator combination. Hydraulic fluids , such as oil. 
even though far less compressible than gas, when compressed 
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in a cylinder compartment act like a spring and introduce a 
second order mass - spring system whose natural frequency 
limits the bandwidth of any hydraulic servo system abruptly. 

At no load the equations of motion of the piston may be 
written as follows: 


= 


= A 


d X 


d t 


4 ft 


d P, 


d t 


ca. aa:> 






ca. ab:> 


p + P^ 
1 a 


ca. acD 


A p, = 
p 1 


m 


d‘^x 


d t 


ca. a:) 


where, = piston displacement 

ft = bulk modulus of the fluid = - 


dP 


dV/V 


= load flow 

= flow rate into chamber 1 

CLs = flow rate out of chamber a 

Pi =• pressure in chamber 1 

= pressure in chamber 2 
c? 

m = mass of bhe piston 
P 

A = area of the piston 
P ^ 

V = total volume of the cylinder 
The first term on the RHS of Equation 2.2a is the rate of 
change of the volume of the two chambers of the hydraulic 
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cylinder, while the second term represents the rate of change 
of volume of the hydraulic fluid due to finite 
compressibility of the hydraulic fluid. 

2.2.2 VALVE DYNAMICS 

Constant supply pressure and sump pressure are assumed. 
The load flow through the valve may be expressed!! 1 3 as a 

function of valve position and the load pressure P^. 

^1 " ^ ^ • ^1 ^ 

Hie above relation may be linearised to yield an 
expression for incremental value ACL for load flow Q 


AQ, 



V 

The r equi r ed 
dif f erentiation 
char acter i sti cs 
curves. 
Defining, 


^1 

Ax + = AP, C2. 43) 

- ^P^ ^ 

partial derivatives are obtained by 
of the equation for the load flow 
or graphically from a plot of load flow 


a 

valve flow gain k = C2. 53) 

V 

valve flow pressure coefficient k = - C3. 6D 

From Equations C2. 43) , C2. 53) and C2. 63) , neglecting the higher 

order terms , we get 

AQ, = kAx - kAP- 
^ g V cl 


C2. 73) 



is negative Tor any valve conriguration v/hich niakes 

the flow pressure coefficient always a positive nuniber . 

The valve flow gain k indicates the increment in load flow 

g 

due to incremental displacement of valve* at a given load 
pressure P^. k^ is always positive for the definition of 
and Q, . The valve pressure coefficient k^ indicates the sum 
of decrease in inflow into chamber 1 from supply* due to 
increase in pressure P^ in chamber 1 and decrease in outflow 
from chamber S to the sump due to decrease in the pressure 
in chamber 2. This also accounts for decrease in flow rate 

1 

and increase in due to leakage across the piston due to 

, the difference in chamber pressuresC 1 3 . 

2.3 LINEAR ANALYSIS OF ELECTRO - HYDRAULIC ACTUATOR CNO LOAD> 
Considering Equations C2. 2Z). C2. 3D and C2.7D we obtain an 
e>q5ression for the transfer function between the displacement 
of the valve and the displacement of the piston of the 
hydraulic actuator 


X CsD 
P 


k / A 

g p 


C2. SD 


X CsD s C s^/ (ii, + 2 l / (j>, s+ i D 

V h h 


4 ^ a! 


wher e c>>^ 


m V 
P 


2 C/ 


k m 
c P 
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For- -the typical actuator considered the hydraulic 

frequency is 6BO. 25 rad/s and C the damping coefficient is 
0.01 2S7 for the nominal values of the parameters of the 

hydraulic actuator. 

The nominal values of all the parameters of the system 

considered for illustration of the analysis as given in 
Appendix B are based on the product information brochure of 
an E H actuator manufacturer. 

Tlie frequency response corresponding to the open loop 
transfer function of the hydraulic actuator is shown in 
Figure 2. G . It may be seen from the figure that at low 

frequencies^ the amplitude ratio has an attenuation at the 

rate of 20 db^^decade; has a sharp resonance at around G50 
rad^^'s and an attenuation of GO db/decade at high frequencies. 

Thie frequency response of the electro hydraulic actuator 
C configuration : DDV + actuator D to input voltage at no 
load» is given in Figure 2.7. The bandwidth of the system is 
SO rad/^s. For an elevon actuator the frequency range of 
interest is usually up to 10 Hz C 65rad/s0. The bandwidth of 
the E H actuator system is therefore adequate for the purpose 
of actuation of the elevens. The corresponding transfer 
function is given in Appendix ACA3D. It is of considerable 
interest to study the changes in the dynamics of the E H 
actuator as some important parameters are varied. 

2. 3. 1 EFFECT OF VALVE COEFFICIENTS 

In the above analysis the valve flow dynamics was 
linearized about a nominal operating condition. C vide Section 
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coefficients may be at variance with the nominal values. It 

therefore seems important to study the effect of the valve 

coefficients. From Equation C2. SD it may be seen that the 

valve flow aain k directly affects the aain of the hydraulic 

- g ' 

actuator. It therefore affects the stability of the closed 

loop E H actuator system. Figure 2. 8a shows the root locus of 

the E H actuator for variation in the value of k . The roots 

g 

move towards the unstable half of the s-plane as k 

g 

increases. Figure 2.8b shows the frequency response of the E 

H actuator for three values Chalf, nominal and double! of k . 

g 

It n\ay be seen from the figure that the resonance frequency 

does not change with the variation in the valve flow gain k^. 

The dampina reduces with the increase of k . 

g 

The valve flow pressure coefficient directly affects the 
damping coefficient. This may be seen from Equation C2. 8! as 
well as from Fiaure 2.9. For half the nominal value of k 

c 

the roots are almost on the imaginary axis indicating a very 

low value of damping. However the increase in damping is not 

substantial as the value of valve flow pressure coefficient 

k is increased to double the nominal value. The resonance 
c 

frequency again remains unaffected by change in k^. 

2.3.2 EFFECT OF BULK MODULUS 

The effect of bulk modulus of the hydraulic fluid on 
E H actuator is studied. The variation of the eigenvalues as 
the bulk modulus {3 is varied from half the nominal value to 
double the nominal value is shown in Figure 2. iOa. It may be 
noted that the system is unstable for low values of bulk 
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modulus. The usual value of bulk modulus for hydraulic fluids 
used in actuators are of the order of S50000 psi . However* 
the nominal value of lOOOOO psi has been arrived at to take 
ii’ito account the flexibility of the cylinder walls and the 
pipe lines which supply the fluid from the sourceC13.It may 
be seen that both the resonance frequency and the damping 
ratio increase with increase in f^. This may also be seen from 
Figure £. 10b Cfrequenc^/ response plots for varying values of 

2.4 COULOMB FRICTION [1,23 

The factors affecting the damping of the motion of the 

hydraulic actuator are viscous friction* intersurface 

friction between the piston and cylinder, which has the 

nature of Coulomb friction, and leakage flow. A difficulty 

here is the non-linear behaviour of Coulomb friction, making 

conventional linear analysis impossible. It has been 

demonstrated by Edwards [63 that a valid model may be assumed 

by setting Coulomb friction equal to zero and introducing Cor 

increasing viscous friction to provide the required daBipi ng 

for a frequency range. A similar model is chosen in the 

present work for linear analysis purposes. 

To justify such a model a digital simulation is carried 

out. Coulomib friction F is simulated as follows CFi cure 2. ii:> 

c 

X 

F = F £ i X X O 

c cmax I • I P 

X ^ 

\ p\ 



F is the maximum value of Coulomb friction between the 

cmax 
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actuator piston and the actuator chamber. 

In physical reality* Coulon±> friction at x = 0 is 

p 

deter itiined by the fact that it will take whatever value is 

necessary to keep the velocity zero. As long as the net 

force* F Charring Coulon±> frictionD acting on a stationary 

piston is - F < F < F , the speed will reniain zero. 

^ cinax ciiaax ^ 

As soon as F exceeds F the piston starts movina. 

cmax ^ 

Figure S. i2a shows a simulated time response to unit step 
input of voltage to the actuator system with and without 
Coulomb friction. A fourth order Runge-Kutta method is used 
to solve the associated differential equations. Viscous 
friction coefficient is now chosen by trial and error method 
to achieve a time response as close to that of the system 
with coulomb friction as possible. The two time responses are 
compared in Figure 2. iSb. The figure shows a very good match. 
The two responses are indistinguishable after 0.03 seconds. 
This appears to justify further analysis of the system using 
equivalent viscous friction for the model. 

Ai^ analysis of the system by Describing Function CD FD 
method throws a different light on how the coulomb friction 
operates. 

In the method of the Describing Function all signals are 
assumed to be harmonic. Assuming the velocity input signal to 
be sinusoidal* the Coulomb friction output function n\ay be 


written as; 
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input : X = X sin a> t 

P P 

4 1 

output : F = F C sin a> t + sin 3 a> t 4- . . . 

^ an cmax 3 

It is important to note that in witing the above 

expression, it has been assumed that the nonlinearity does 

not generate sub harmonics. As the Coulomb friction function 

is an odd function of x the even harmonics terms do not 

P 

appear in the output. The absence of a phase shift is due to 
the memory-less nature of the Couiojnb friction function 

[Q] . 


In the DF - method one takes the first harmonic only, 
ignoring the higher harmonics. So, the DF - method gives a 
linearized relation between the input and the output through 
an equivalent transfer function. 


F p 


F Cc*>:> 
c 


X CoD 

P 


4 


n 


F 4 F 

cmax cmax 


X fr a> X 

P P 


ca. s:) 


Putting the above linearized transfer function in the loop 
and taking the root locus by varying K- it is clearly seen 
from Figure 2.13 that as Kp varies from O to co, the damping 
factor of the dominant root increases. Fi'om the above one may 
infer the following: 

iD Coulomb friction contributes positively to the damping of 
a hydraulic actuator. 


i i 5 Let 

F be of a 

cmax 

constant value. 

I n that case as 

X 

P 

varies 

from O to 00 , 

Kp varies from 

0 

8 

O. Thus 

at 

a gi ven 

frequency, as 

di spl acement of 

the 

pi ston of 

a 


hydraulic actuator reduces, the coulomb friction contributes 
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more tov/ards the damping of the system. 

The above inferences made from the describing function 

study of the Coulon\b friction function are proven correct by 

actual simulation of the electro hydraulic actuator system in 

which the Couion\b friction function is directly incorpor ated. 

Figure £.14 shows the time response of the actuator to step 

inputs of varying magnitudes viz. 1 vdc* 5 vdc and 10 vdc. 

The nominal level of Couion±) friction F is SO lbs. 

cmax 

Compared to the response of the system with no Couiontib 

friction* which is lowly damped* high frequency oscillatory* 
as sho'wn in Figure £. i£. Coulomb friction uniforinly damps 
these oscillations. 

Furthermore* it may be observed from Figure £.14 that the 
damping due to Coulomb friction reduces as the magnitude of 
the response C displacement 2) increases. 

In addition* a study is made of the effect of the level 
of Coulon:\b friction on the response of E H actuator which 
is shown in Figure £.15. It may be seen that as Coulomb 

friction is increased from 80 lbs. to 8000 ibs, the damping 

increases* resulting in increased settling time, 

£.5 MOUNTING STIFFNESS 

The actuator so far , has been assumed to be rigidly 
mounted on the airframe, A flexible mounting is now 
considered. The pressure developed in the cylinder chambers 
now not only accelerates the piston but also moves the 

cylinder in the opposite direction. 
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The Equation C2. 2aD now gets n^odified as follows: 

^ " ^a ^ V d P, 

- + L_ c2.io:> 

^ 4/5 d t 

where is the displacement of the cylinder and is given by 


“ . 

where kj^ is the mounting stitfn^gg expressed as an equivalent 
spring stiflness. It should be n^ej^tioned here that the value 
of cylinder mass for the actuato^ considered for illustration 
is not available. However, simf^t^ calculations show that at 
the operating frequencies the inertial force associated with 
possible values of this mass insignificant compared 

to the force due to the mountij^g stiffness at its nominal 
value for the system considered. therefore the mass of 

cylinder may be ignored. 

The transfer function betweej, the voltage input to the 
DDV and the actuator displac^^^^t, and the corresponding 
eigenvalues are given in Appendtj^ A C A4 5 

Figure 2.16a shows the root locus for the variation of 
the mounting stiffness from nominal value to 

double the nominal value. The nearer to the 

imaginary axis show an increase natural frequency as well 
as an increase in damping wh^^ stiffness is increased 

from O.l9<nominal valued to 1.5*Knoininal vaiueD. Beyond that 
the damping decreases and the hatural frequency increases. On 
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the other hand the roots farther from the imaginary axis show 
a different trend. The natural frequency decreases as the 
stiffness is increased. The damping first increases, reaches 
a peak for a stiffness value of 0.7 times the nomnal value 
and then decreases. 

Figure S. 16b compares the frequenc2/ response of a 
flexibly mounted E H actuator at the nominal value of 
stiffness with that for rigidly mounted actuator. The 
bandwidth of the E H actuator with flexible mounting is 140 
rad.xs Ccompared to SO rad^s for the rigidly mounted 
actuatorO. However, at higher frequencies the attenuation is 
greater for the E H actuator with flexible mounting. 

The step response of the E H actuator is shown in Figure 
S. 16c for a finite Cnomlnal value!) and infinite value of the 
mounting stiffness. It may be seen from the figure that the 
damping has reduced and there is an overshoot of about 11% 
when the mounting stiffness is finite. The rise time has 
however reduced as is also indicated by Figure 2.16b. 

A frequency response study is carried out for variation 
in mounting stiffness to understand this phenomenon, and 
shown in Figure 2. 16d. The bandwidth is seen to be increasing 
as the mounting stiffness is reduced from infinity. The 
maximum is reached at the nominal value of the mounting 
stiffness. Any further reduction in the mounting stiffness 
results in a reduced bandwidth. The explanation may be as 


follows: 
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eigenvalues of E H actuator Crigid mount! ngD : 

-109.87 ± 641.58 i Coi = 650.86; C = 0.168S:> 

-119.36 ± 14.58 i C(/" = 120.84; C = 0. 9926D 

eigenvalues of E H actuator Cflexible mounting!? ; 

Cfor nominal values of the mounting stiffness!) 

-66.58 ± 122.66 i Coi = 135.56; C = O. 477D 

-154.98 ± 80.09 i Co/" = 174.45; f = 0, 8SSD 

n ^ 

From the above values it may be seen that the value of 
the undamped natural frequency of the lov/er frequency mode 
the E H actuator on flexible mounting is greater than the 
cor responding value for rigidly mounted E H actuator. 
Similarly the damping factor for the E H actuator with 
flexible mounting is lower than that of the rigidly mounted E 
H actuator. Both these factors contribute to a larger 
bandwidth. 

As the stiffness is further decreased to half the nominal 
value, the eigenvalues of the system are as follows: 

-29.07 ± 92.19 i Co^ = 87.18; C = O. 3334:> 

n ^ ^ 

-192.11 ± 63.69 i Ca> = 202.39; C = 0.9492> 

n 

It may be seen that the lower end natural frequency value 
is now less than that for the nominal value of the mounting 
stiffness. Besides, the damping ratio is more. Hence a 
smaller bandwidth fiiay be expected. 

2-6 ACTUATOR - LOAD DYNAMICS 161 

Figure 2.17 shows schematically the actuator - eleven 
servo system. The elements included for the analysis of this 


system apart from the DDV, hydraulic actuator- 


and the 
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mounting stiffness are: 
iD Inertia of the eleven CI3 
iiD Horn stiffness 
iii^ Aerodynamic Hinge Moment 

It is assumed that the elevens are rigid. The flexibility 
of the linkages C piston rod and eleven horn D is represented 
by a linear spring at the horn coaxial to the piston rod. 

Since the n:^ass of the cylinder is not accounted for, the 
whole system is represented as a two niass system, whereas, 
strictly it is a three mass system. Figure 2.18 shows the 
siejn convention of the displacements, linear and rotarv, of 
the actuator - elevon system. 

The equations of motion are as follows: 

''h ="p ■ ^ ^ ‘ ^ '‘<5 * '‘o “ 

whisro 1 == horn leng'th 

6 = elevon deflection 
OL = angle of attack CAOAD 

k- = hinae moment per unit elevon deflection 
a 

k = hinae moment per unit AOA 
a 

Hie angle of attack inputs may be considered as 
disturbances for the present and the transfer fui/|Otion 
between the piston displacement and the elevon deflilition 
follows from Equation C2. 13D 

6Cs> k^ 1 

I s^ + C K. 1^ + k c D 
h o 


X CsD 

P 


C2. 14j| 



Equa_iLon 2.12 is shown in Figure 2.1 Ccategorized as load 
dyrianii cs 3 . 

The o-ve-rail transfer function of the actuator *- eleven servo 

systeeii:! c^bhalned from Figure 2.1 is given in Appendix A CA5D. 

The ei genvalues are : 

-100.72 ± 8S4.41 i 
—234, 21 
-2S. 36 

■+10. 44 ± 28. 63 i 

The load resonance poles have positive real part. The 
SYslea\ is therefore unstable for the nominal position 
feedback: ^ain. It is instructive to study the root locus of 
the E H actuator for varying liiass. Figure 2.19 shows the 
doim rant roots moving towards the right half of the s - plane 
with Inoreasing mass C the mass considered is the mass of the 
pistoni ai-^d the equivalent mass of the eleven due to its 
mofnent of inertia. The two masses are assumed to ba rigidly 
connected to each otherD. The roots become unstable for 
an eguiivadent mass greater than 3. O lbs. 

The instability may be due to the over cor recti on of the 
posi tLor feedback. Tine feedback gain may have been determined 
for a no “ load case and hence the instability appears due to 
eievoin inertia. A possible solution is the damping of the 
load poJLes by suitable means. 

2. 6- 1 DjARPING by dynamic PRESSURE FEEDBACK: 

Quite a few methods exists for damping of hydraulic 
servosC2] . Prominent among them being: 

a) Dampi ng through increased leakage C disadvantage : 
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hi gh pumpi ng r equi r ementsD 

bl) Damping through pressure feedback C disadvantage 
increased load sensitivity D 

cl? Damping through acceleration feedback C di sadvantage : 

increased hardware complexity D 
cD Damping through Dyi’iamic pressure feedback CDPFD : 

Dynamic C or derivative 5 pressure feedback is usually 
preferred because this does not change the load 
stiff ness. 

It is represented in the transfer function form associated 
with what is known as a washout circuit. 


K 


m 


k 


P 


T S 

P 


T S + 1 
P 


C2.is:> 


At low frequencies = 0 and the load stiffness therefore 

remains unchanged. At higher frequencies. will 

provide adequate damping if oj t » l.The block diagraai of the 
actuator - elevon servo system with Dynamic Pressure Feedback 
is given in Figure 2.20. 

T is determined such that it satisfies the condition u>t»1 

Thus T = 3 is taken Cdetermined by trial and errors. Root 

locus is then drawn, as shown in Figure 2.21, to deter nine 

the value of k . k is varied from 0.0002 to 0.001. Best 
P P 

results are obtained for k = 0.0006. The transfer function 

P 

of the actuator - elevon servo system with dynaaxic pressure 
feedback is given in Appendix ACA6D. Figure 2.22 shows the 
unit step response of the actuator - elevon servo system for 
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T = 3 and k = O. OOOS. I*t mav be seen from the Fiaure that 

p 

the transient response is oscillatory with an overshoot of 
about 66%, The settling time is about two seconds. 

As may be seen from Appendix ACA6D, the damping ratio of 
the dominant load resonance poles is 0.113. Evidently the 
dynamic characteristics of the elevon actuation system are 
poor although those of the E H actuator under no load 
conditions were satisfactory, as shown in Figure 2. l£b. It 
ftiay be noted that the above character! sties of the elevon 
actuation system were the best Cin terms of rise time and 
settling timeZ? that could be achieved by adjustment of the 
two parameters and of the dynamic pressure feedback 

loop, the rest of the parameters of the system considered for 
illustration being preserved. The E H actuator appears to 
have a tnismatch with the impedance char-acteristies of the 
eleven system. 

While one may attempt at changing the structure of the 
system by introducing more feedback loops, or changing the 
controller of the dynamic pressure feedback loop, which has 
presently the form of a washout circuit having two 
parameters, it may be worthwhile to see how the system 
performs, when placed as it is in the longitudinal flight 
control system of the aircraft. 

2.6.2 EFFECT OF VARIATION OF HORN STIFFNESS 

A study is made to show the effect of variation of the 
horn stiffness on the elevon actuation system. Figure 2.23a 
shows the root locus of the dominant poles when the horn 
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stiffness is vat'ied from 10% of the nominal value to doable 
the nominal value. It may be seen from the Figure that the 
sensitivity of the roots to variation in horn stiffness is 
very little above stiffness value of 0.7 times the nonuLnal 
value. The S 3 /stem natural frequency and the damping factor 
reduce for sii\alier values of the horn stiffness. The same 
inferences may be dravm from the frequency response curves 
which are showi'i in Figure 2.23b. 
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CHAPTER 3 IMPEDANCE OF THE ELEVON ACTUATION SYSTEM 

3.1 INTRODUCTION 

The restraint to displacement of a control surface which 
is provided by its control system may be termed as its 
’ impedance^ C 8] . The stiffness of the actuator loop must 
consider all elements including stiffness of the horn on the 
control surface, bending of bolts, internal stiffness of 
bearings, free pia^*^ and friction in bolted connections, 
actuator structural stiffness, actuator hydraulic stiffness 
and spring rate of the tie down structure etc. [43. 

3.2 IMPORTANCE OF IMPEDANCE CHARACTERISTICS STUDIES 

Impedance characteristics of the eleven actuation system 
are usually determined to carry out aeroser voelastici ty 
studies. Though aeroser voelastici ty is beyond the scope of 
the present work, it may be worthwhile to determine the 
impedance characteristics of the elevon actuatioi^i system as 
an input to further studies on flutter control of the wing, 
the control surfaces etc. 

Table 3.1 reveals the importance of impedance studies 
[43. It summarizes potential solutions for the different 
types of conbroi surface flutter C different from wing flutter 
and characterized by higher frequency of 40-S0 Hz> with 
comments regarding relative effectiveness of each solution. 
Manipulation of the control system stiffness appears to be 
the most efficient means of flutter supression. 
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TABLE 3.1 

SOLUTIONS TO WING FLUTTER CONTROL 




CANDIDATE 

SOLUTIONS 


FLUTTER 

TYPE 

MASS 

BALANCE 

HYDRAULIC 

DAMPER 

RESTRAINT VTA 
STIFFNESS 

HINGE LINE 
LOCATION 

Type 1 

classical 

coupling 

i 

very 

err ecti ve 

use is 

rare 

ALMOST 

ALWAYS 

THE MOST 

bai anced 
surfaces 
present less 
pr obi em 

Type £ 

al 1 mo vabl e 

surfaces 

1 i mi ted 

1 

t success 

at ti mes 

1 

use is 

rare 

1 

1 EFFICIENT 
SOLUTION 

IN TERMS 

j 

ef f ect 

difficult to 
gener ai i ze 

] 

Type 3 
transonic/ ; 
supersonic | 

j 

not 

ef f ecti ve 

very 

effect! ve 

1 

OF 

LBS^'KNOT 

IMPROVEMENT 

not a 
major 
variable 
















The dynamic effect of actuator impedance is illustrated 
in Figure 3.1. Here Z C *ia> is the mechanical impedance 

of the actuator on a flexible mounting, expressed as 
equivalent hinge moment per unit eleven rotation. is the 

steady aerodynamic hinge moment appropriate to the flight 
condition ux'^der consideration. If feed back control is not 
employed C i.e. the open loop case D to suppress flutter we 
have [7] 

2 = Q. 

where is the amplitude of elevator deflection 
and expCiootD is the total dynamic eleven hinge 
moment viewed as a loading action. 

Thus from Figure 3.1 we get 

q = T Ca d + 

^ ^ z 

where T is the transfer function of the 

actuator elevon servo system under steady 
loading and d is the demand signal. 

The above two cases illustrate the need for determining 
the impedance of the actuator elevon servo system. Linearized 
results from actuator impedance tests are sometimes used. But 
the impedance may also be estimated by mathematical 

modelling of the actuator. 

3.3 ACTUATOR HYDRAULIC STIFFNESS 

Impedance of a mechanical system depends on clearly 
defined physical entities viz. spring stiffness, dampers. 
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masses and their coinbi nations. Hydraulic actuator stiffness 
however depends on factors like piston area, bulk modulus, 
servo action provided by feedbacks etc. to provide the 
necessary impedance. A simplified hydraulic actuator Csingle 
cylinder; rigid mounting; Figure 2. 5D is analysed to 
understand the actuator hydraulic stiffness. 

Let the direct drive valve be represented by a unit 


gain. 


Then the equations of motion are: 

o V cl 


= k_ x_ - k. P, 


d P, 


O = A k + 
^ p P 


4/5 d t 


i = V - k X - k 
V o X p p 


d P, 


d t 


C3. lal 


C3. lb:? 


C3. ic!) 


A = F C3.1d5 

where k = position feedback gain 

k = dYi'iaraic pressure feedback aain 
P ■ 

V = conmvand voltage 
o 

other terms are the same as used in Chapter 2. 


Pqj-- V =0, from Equation 3.1a - 3. Id we get, 
o 

FCsZ) CApS + <^4/5Ap/ VsD 

Cl + C4/Sk / Vsl? + C4/5k k s^'Vs:>> 

c g p 


X Cs> 

P 


C3.2:) 
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3.3.1 CASE Ai No feedback to the valve Ck^ 
Equation 3.2 reduces to 


0; k = OD 

X 


FC sU 



X Cs5 
P 



k 


C 3. 33) 


If* we assume for the present that the hydraulic oil to 
be incompressible C.3=ai>, Equation 3.3 reduces to 


FCs:? 


SXp^So 



C3. 43) 


Physically Equation 3.4 n)ay be explained as follows; 

The nxition of the actuator piston induces load flow which is 
by the finite opening of the valves and possible 
leakage across the piston, as reflected in the inverse of k^. 
tbe valve flow pressure coefficient at Equation 2. © . 

Now let the hydraulic oil be considered to be 

cojnpr essi bl e. 

£jQf"ining effective oil spring k^£S3 a^s 


4 /? 


k = 
o 


P 


v/e obtain from Equation 3.3 

FCsy _ 


k s 
o 


X Cs3) 
P 


s + k k- / 

o c p 


C3. S> 



31 


For high values of input frequency this reduces to 


FCs=jt>>D 


X Cs=JoiD 
P 


From the above equation it may be seen that at high 
frequencies the displacement of the actuator piston is 
essentially resisted by the effective oil spring. An increase 
in the bulk modulus increases the stiffness of the hydraulic 
actuator. For low frequencies the expression for impedance 
reduces to Equation 3.4 which indicates resistance to motion, 
not to displacement. 

Alternatively equation 3.5 may be viewed as follows: 

If flow pressure coefficient k is very small such that the 

C " 

2 

term k k /A mav be neglected Ci.e. the load flow is 
o c p " 

negligibie!> the resistance to motion of the piston is 
provided by the effective oil spring. For high values of k^, 
the load flow alone determines the amount of resistance which 
is provided by the actuator. 

3.3.2 CASE B: k = O; k O ; ;? = oo 

p X 

In this case the effect of position feedback is 
investigated. Equation 3. 2 reduces to 


FCs:> 


X CsD 
P 


k k A 
g X p 


A”^ s 


C3. ey 


The second term has been already investigated. For lower 


values of frequency the impedance value tends to be 
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_ k k A 

FCs=ja>3 _ g X p 

X Cs=ja>D k 

p c 


Reference [83 defines Ck A y'k D as effective actuator 

g p c 

servo stiffness. Tl-^us* for lower values of frequency/ 


Actuator hydraulic stiffness = position gain feedback * 

effective actuator servo stiffness 


Physically, what happens is as follows: 

Without loss of generality it may be assumed that 
initially the valves are closed. A change in the piston 
position causes the servovalve to move in the opposite 
direction due to the position feedback and the consequent 
reversal of flow of hydraulic oil resists the displacement. 
3.3-3 CASE C: k O; k = 0; k =0 

In this case the effect of dynamic pressure feedback is 
investigated. Equation 3. 2 reduces to 


A 

s FCsD ^ p 

s X Cs::) C k k V + D 

P 9 P 


C3. 6D 


Equation 3.6 may be interpreted as follows: 

A motion of the piston in positive direction CFigure 2.53) 
results in outflow of hydraulic oil from chamber 2 along with 
the compression of the oil in the chamber and associated 
rise in pressure. In chamber 1 there is an inflow along with 
drop in the chamber pressure and expansion of the oil. Thus a 
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constant velocity input results in a rate of change of load 
pressure “^ 2 ^ which is fed back to the valve through 

the dynanvLc pressure feedback which results in a change in 
valve poslt-ion so that the outflow from chamber £ and the 
inflow to chamber 1 reduces thus providing resistance to the 
the motion of the piston. 

3.4 DETERMINATION OF IMPEDANCE OF THE ELEVON ACTUATION 
SYSTEM 

Initially the n\ass is not considered. This enables us to 
determine the actuator hydraulic stiffness. From Figure 
2. 1 we get , 


eCs:) = V - kda^klaJ<X CsD - k - 
o p p 


T S 


T S+1 

P 


CP, , CsD+P_CsDD C3. 75 
11 12 


eCs5 fCT 
I T s + 1 

eCsD KT 
t r s + 1 


F, Cs5 - Ap A 


s A X Cs5 
P P 


s A X Cs5 
P P 


- s A 


Pi, J 


Vs + 


F, Cs5 
1 


C3. 85 


F Cs5 4,3 A 

s A — ^ ^ = F..,Cs5 

p . I ,, ...... 2 


k. 


Vs + 4/3 k 


b ' c 

C3. 95 

wher e KT = ka»<kl ^rvs^k 


P, , Cs5 A = F^ Cs5 

11 p 1 

P Cs5 A = F_Cs5 C3.105 

X ^ |3 cS 

F. Cs5 = F^Cs5 


NOTE: 15 In the above equations eCsD is the Laplace transform 
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of the error signal to the system. 

The subscript 1 & 2 refer to the tv/o cylinders of the 
dual “ tandem- actuator . 

3D It is assumed that the force generated by the two 
actuators is identical. 

From Equation C3. 8D. C3.QD & C3. lOD we get, taking V = O 


FCsD 

X CsD 
P 


[sCTS+lDe/^A*^ + 8KT kda kla ftA 3Ct s + ID 
p p p 

Cts+1>Ct s+i:)CVs+4/?k sCts+15Ct s+1D + 

P ^ ^ b ^ 


C3. 113 


CKT 8/? k T 3s 
P P 


Now considering the mass and equivalent damping we get 

FCsD = FCs5 + Cm + B sD X Cs!) 

P V p 


FCs!) FCs:) 2^_ . 

= +Cm s +B sD 

X CsD X CsZ) ^ 

P P 


C3. 12:) 


Actuator impedance C D and the horn stiffness may now be 


considered as two springs in series. 


Hence, the equivalent Impedance is 



C3. 133 


Impedance is usually expressed in terms of required torque 
per unit deflection of the control surface 


Therefore, 
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T _ , 2 

— i ^ 


eq 


C3. lO 


The eleven may now be considered to be restrained by a 
torsional spring. The torque required to deflect the eleven 
has to overcome the torsional force due to this spring, the 
hinge moment for a given flight condition and the Inertial 
force of the elevon. 

Hence, 


T = T + I 6 + H C3. is:) 

s 

wher e, H = k 

s o( c5 

For a given flight condition considering the hinge moment due 
to elevon deflection only we get the following closed form 
sol ution: 


TCsD TCsD 

= -1- c I s*^ + k- :> 

dCsD dCsZ) 


= l^Z ->• CIs^ + k<,> 
eq o 


C3. lb!) 


Figure 3.2 shows the impedance characteristics of the 

actuator elevon servo system versus frequency for the 

following flight conditions: 

i5 altitude = 5 km ; Mach No = 0.3 
k^ = 47. 58 N m / deg 

iiD altitude = sea level ; Mach No = 1.2 
k^ = 4627. 85 N m / deg 

iii> Aircraft at rest k\. = O 

o 

The following observations may be made from Figure 3*2. 

a> The dip In the impedance curves is at the load 
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resonance frequency Ci.e. 15.9 rad/s for flight condition Ci:> 

and 32.6 rad/s for flight condition Cii2). The impedance 

curve for aircraft at rest overlaps the curve for flight 

condition CiD. This is because of the small difference in the 

aerodynamic elevon stiffness k- values. 

o 

b2) For low values of frequency, an increase in 
aerodynamic elevon stiffness predictably causes an Increase 
in impedance of the elevon actuation servo. However for 
higher values of frequency the trend reverses. 

c> For ->■ 00 , the impedance of the elevon actuation 
servo system does not change much with the change in flight 
condition. This is due to the fact that at higher values of 
frequency the coefficients of s and higher powers of s in 
Equation 3.16 predominates and impedance of the elevon 
actuation system is thus no longer sensitive to the variation 
in k' . 

Figure 3. 3 shows the phase plot for the impedance of the 
elevon actuation system. It may be seen from the figure that 
for lower values of frequency the elevon deflection is almost 
in phase with the external loading. As the frequency 
increases the phase difference tends towards 180*. However, 
for intermediate frequency values C16 to 32 rad/s 
approximate! y> the phase difference between the external 
loading and the elevon deflection is greater for low values 
of aerodyriamlc stiffness compared to that for higher values. 
The relative phases of applied torque and elevon deflection 
have a strong bearing on the flutter characteristics of the 
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control surface. 

3*5 RELATION BETWEEN IMPEDANCE AND THE TRANSFER FUNCTION OF 
THE ELEVON ACTUATION SERVO 

To determine the relation between the impedance and the 
transTer function C= control surface deflection / cofiuiiand 
signalD of the eleven actuation system, it is necessary to 
consider the inverse of impedance, i.e. compliance. 

Compliance may also be defined as the transfer function 
between the control surface deflection and the external load 
C compliance = control surface displacement external 

torqueZ). Figure 3.4a Cderived from Figure 2. 1Z> shows the 
relation between the transfer function and the compliance of 
the eleven actuation servo. If the conunand voltage is 
considered to be zero. Figure 3.4a gives the expression for 
compliance. If the external load is considered zero, the 
figure gives the expression for transfer function between the 
command signal and the eleven deflection. 

Thus , 

ea G3 

compliance = £ = 

FCsD 1 + G1 Ga Ha + Ga G3 HI 

X Cs> G1 Ga G3 

transfer function = ^ = 

V Cs> 1 *4. G1 Ga Ha 4- Ga G3 HI 
o 

The simplified block diagram representation of the 
eleven actuation system is shown in Figure 3. 4b. It may be 
seen from the Figure that 

Ga* 

compl i a nee = ; r 

14- GiGa Hi 
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Gi Ga 

transfer function = ; , 

1 4* G1G2 HI 

Figure 3. 5 compares the frequency response of the 
transfer function and the compliance for flight condition 

CiD. The two curves show similar behaviour. The response is 
uniform up to 10 rad-^'s, with a resonance peak at around 16 
rad/3. 

Figure 3.6 shows the difference between the transfer 

function and the compliance in decibels. The transfer 

function Gi4<l is also plotted. Gl may also be seen as the 

transfer function between the force in the piston rod and the 

conm\and voltage input as the piston position is held 

constant. Opening the valve through input command voltage 

results in changes in and consequently in the force 

developed in the piston rod of the actuator. The position 

feedback to the DDV will be ineffective when the actuator 

position is held constant. eventually reaches the value 

P , - P , for a constant command voltage . 

supply/ sump 

Compliance indicates the effect of disturbances* 
including the aerodynamic forces due to gusts, wake or 
flutter* vibration of the aircraft and control surface* on 
the response of the control surface, as the command signal is 
given to the actuator system. Thus compliance viewed as the 
sensitivity of the control surface deflection to forces 
considered as noise/'^disturbances* should be as sii^aii as 
possible* i . e. impedance should be high. 

It nm.y be noted from Figure 3.4b that the forward path 
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wit.h .=»plidnc« i= sa-. while theh ee=eelel..d wlhh 
l.eeefer fenehlee between end d 1= Glsa'. SI end GS 

.ne eleeely neleted In the eenee thet ohenglng pereeetene ef 
G2* t.o eomplianea is not possible in most casos 

without aimultainoouBly changing G1 . G1 and G2 ai'* 

i rici^p«!?rnd«?rit/ pHycaicAl 45rrit#i‘Li<ft’®- 

Figure 3.4c shows the simplified block diagram 
representation of the impedance of the eleven actuation 
system. From the figure it may be seen that a deflection of 
the control surface is resisted by the control system through 

the following two forward paths: 

a5 GlHl * ; which represents the transfer function 

b^^tween the resistance offered by the actuator due to its 
servo action through the elevon position feedback and the 


tfbX OtWaOvint iOtO 'Lr X ArB . 

bD 1/GS^ ; which repri&sen-Ls the tramsrer function 

between the resistance offered by the inertia Cmass of piston 
and control surf ace!) > Coulomb frietion» aerodynamic and horn 
stiffness etc. to the elevon deflection and the elevon 
deflection. This includes the resistance offered to motion by 
the hydraulic fluid. 

The impedance may be thus represented as 


resisting load TCs3 


1 + GiGS^Hl" 




elevon deflection 


<5Ca:> 


V =o 
o 
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CHAPTER 4 ACTUATOR AIRCRAFT INTEGRATION 


4.1 INTRODUCTION 

Electro hydra.ullc servo mechanisms are used to provide 
power boost, stability augmentation control inputs, or 
primary control of the airplane. The present generation 
airplanes are control configured vehicles wherein the static 
stability requirement is relaxed at the design stage in order 
to increase maneuverability. ' In such airplanes the actuator 
forms an integral part of the overall control law scheme. 

4.E AIRPLANE MODEL [ 3, lO 3 

The airplane considered here for illustration is a 
fighter aircraft of relaxed static stability. The primary 
control loop of the longitudinal control system consists of 
pitch rate and normal acceleration as feedback variables to 
the longitudinal control surface actuators as shown 
schematically in Figure 4.1 . 

Since the system is a Single Input Multi Output C SIMO D 
system, the airplane is represented in State Space form : 

X = A X + B U 
Y = C X + D U 

Where X = state of the system 
A = State Matrix 
U - Input vector 
Y = Output vector 
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The linearised small perturbation equations of longitudinal 
motion of the airplane are: 


m 

[ 

u 4* q w 

r 

V 

3 

= d 



C4. 13 

m 

[ 

V/ 4- P V 

“ q 

u 

3 

= d 

^F 


C4. 23 

q 

I 

4- D r C 

I 


I 

3 

- r^ I 

+ p^ I = d M 

C4. 33 


3" 

X 


z 

xz 

^ xz 



X„ = X„C ot, q. 6 y 
F r e 

Z„ = Z_C a, q, <5 > 

Jh r & 

M = M C a. q, 6 > 

e 

For simplicity, ntaking the short period approximation of the 
dynamics of the airplane, we set u=0 and ignore the equation 
for X„. 

r 

Equations 4.1 ~ 4* 3 reduce to 


C4. 43 


w 


u 


= a = 


m 


a z a z a z 

a 4* q 4- 


a a 


a q 


a s 


C4.Sa3 


( a M a u a H 

a 4- q 4- 6 


a a 


a q 


d 6 


C4. Sb> 
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1 
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we have 


r « I r z 

j I '1 ] [», 


"4 . 


C 4 . G3 


Motina that 


normal accn. at station distant L aft of 


ai r pi ane 


uC o(-q5 + qL. 


c. g. of the 


C4. 7]) 


From Equations C4.5:» and C4. ?:> we get 

'*■ ^ + m l D q + C uZ + M. L > 6 

s H qs o 6 se 

e e 

C4.83 

Using the airplane data available for the flight condition 
Altitude = S km ; Mach No = O. 3 we may write the equations of 

motion Cthe state equation^ and the equation for output 
variables as 


1.701 -0.6174 


-O. 3710 


-7. 1030 


0*0 




The open loop eigenvalues for the svst 


-1 . 9396 


em are: 


O. 6789 
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The aii'piai'ie is thus inherently unstable. The primary control 
loop consists of pitch rate and normal acceleration as 
feedback variables fed to the elevon actuator through a 
filter network and a gain. .The filter constants are: 
i:> Pitch rate feedback filter: 

qfil = 
q gain 

ill) normal acceleration feedback filter; 

1 

a fil = 

^s O. 3S s + 1 

a^ gain = GA = 0.014 
s 

The resultant closed loop eigenvalues are : 

-1 . 9238 ± 1 . 4999 i 
-3. 6044 ± 2. 4465 i 

which are stable. 

4.3 ACTUATOR IN THE LOOP 

The actuator model, as developed in Chapter 2, is put in 
the priniarv control loop as shown in Figure 4.1 using the 
appropriate hinge moment derivatives for the given flight 


O. 31 s + 2. 25 


O. 25 s + 1.0 
= GB = O. 3721 


condi tion. 
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The eigenvalues with the actuator in the loop are: 

-100.77 ± 85S. 28 1 

-118. 99 ± 102. 73 i 

-1.76 ± 15. 90 i 

-1 . 94 

-3. 15 

-O. 24 

-3. 79 

•+0. 57 


The eigenvalue shows jwarginal instability. Tlius it may 
be seen that the system needs some tuning. A tuning gain is 
provided in the feedback loop so that a good response is 
provided. A root locus is plotted for the dominant poles for 
the tuning gain varying from 1 to 15 in Figure 4.2 . The 

figure shows that a tuning gain above 12 may give desirable 
response. Step response of the airplane is shown in Figure 
4.3a Cnornval acceleration vs timel> and 4.3b Cpitch rate vs 
timeD to a unit conmiand voltage for various values of tuning 
gain. It may be seen from the figures that the airplane 
response with a tuning gain of 14 is the smoothest. A tuning 
gain above 14 gives a wavy response while a tuning gain less 
than 14 is comparatively sluggish. Thus a tuning gain of 14 
is chosen. 

The eigenvalues for the actuator - aircraft system with 

realistic actuator with a tuning gain of 14 are : 

-100.77 ± 855.28 i 
-118.98 ± 102.73 1 
-O. 93 ± 15. 57 1 

-3. 57 ±1.62 i 

-1.37 ±1.42 i 



45 


The step response of the airplane with the realistic 
actuator with a tuning gain of 14 is now compared with the 


response 

of the airplane 

wi th an i deal 

actuator 

in Figure 

4. 4a and 

4.4b .It may 

be seen from 

Fi gure 

4. 4a 

that 

air pi ane 

r esponse wi th 

the realistic 

actuator 

model 

is 

si uggish 

cojnpared to the 

airplane response with 

the ideal 


actuator. This may also be seen from the fact that the short 
period frequency of the aircraft with realistic actuator 
model is lower, 3.93 rad/s compared to the cor responding 
frequency of 4,36 rad/s for the aircraft with ideal actuator. 
The aircraft with realistic actuator has a higher short 
period damping ratio of 0.9106 compared to O. 8374 for the 
aircraft with the ideal actuator. 

The overshoot of the normal acceleration has increased 
to 5.35 % from 3% for the ideal case and the rise time is 
higher, 1.0733 seconds compared to 0.875 seconds for the 
ideal case. 

From Figure 4.4b it may be seen that the pitch rate 
response for the realistic model is less smooth, though the 
overshoot has been reduced to 73.93% compared to 91.91% for 
the Ideal case. The rise time is nearly the same - 0.34 
seconds compared to 0.33 seconds for the ideal case. 

The settling time for both the responses Ci.e. normal 
acceleration response as well as the pitch rate response^ is 
seen to have increased for the realistic actuator from 3 


seconds Cfor the ideal case]) to 4 seconds. 
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sl'iitilly is now car'r±eci out L-o d eU er- iiil_ ne if a bet-ter 

per'f cji' Dajuncze raay be obtaii’ied by havi-ftg am eLevon posit-ion 

f'e« 3 d lia. ck ; ili’ir-st-eacl oi' actuiabor' posi ti. ori fe-edBba«.clc:. Figure 4. Sa 

shows Ihoe rcoot- locus or bhe ai r pi a.ne Comity =sh -or t period root-s 

shown) for variation In el evor» feecdbfflc3c gain from the 

raofid naul r-ai ue to S. S times t-he noiiidnal -vaHu-e of position 

Teedliaclc: gadin. It may be seen ff'oni th»e fi_gt_ires -th.at for lower 

vaiues of” eleven feedback ctain the s bo-rU period damped 

freqvie-nccy irs higher and tbe dainpi ng is i es^s. However as the 

jTeedbaLcIc: -gau-in value increases the daajsowd fmecqtxency reduces 

an«d tHe dlaneiping Increases, it is fotjnd liha't for an elevon 

Feed ba-clc: ga«in of 3. S tln^s t-he noiiilnal -^aJ-u e of position 

Feedbau-clc; ^ t-he short period root-s a»ree *,he nearest to 

t-hskt w'iL.h ideal actuator . 

rhe esl»genvai ues of the airplane with ■'ft.al istic actuator 

model fenr l»ie optinium elevon feedback: g-al n (=3. 5»ek^D are: 

-lOO. 79 ± ess. 75 i 
-118. 36 ± 96.34 1 

-1.04 ± 12.73 1 

-3. 6S ± 2. 291 1 

-1.7S ± i. 49 1 

XThe li_i»e response of the aJLr pi arae w-iUh the realistic 

ac-tirat-om iiicodel with optimal eleven . posdt_lo n feedback gain 
cojiipai—ecri ¥_tth that with the nonnixnal. \r-al uess of actuator 
posi tLom rfe-edback and the ideal actuator L,s sfcown in Figure 

-4. Sh» an«d 4 .- 50 . It may be seen from line f~lgur-'es that there 

seeir»s t_o h=)e an improvement in the mesp^nrse with optimal 
el ev'oi'a zl'fr-edBback compared to that vwlth t-has jro_jiii nal position 
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feedback. Compared to the response with the ideal actuator 
the eleven feedback response is wavy. However, the settling 
time is the same. 

It may be erroi-^eous to claim in general that eleven 
jTeediback is better than the position ieedhacki based on the 

above r es ul ts v/hi ch cempar e the ei gen val ues for opti mal 

eleven position feedback with those for nomi nal » non optimal 
value of actuator position feedback. A study is therefore 

made to find out the optimum value of actuator position 

feedback gain k . The optimum k is found to be 0. nominal 
valued. The eigenvalues of the airplane with the realistic 
actuator model with optimum actuator position feedback gain 
are: 

-100.77 ± 955.43 1 
-118.41 ± 101.17 1 
-1.0 ± 12. 5S 1 

-3.69 ± 2. 28 i 

-1.75 ± 1 . 49 i 

Comparing the above eigenvalues with that for optimal elevon 
position feedback shows that the values are very close. The 
response is therefore expected to be the same as there is no 
difference in the short period roots of the airplane Ci.e. 
-3.69 ± 2. 28 1>. This jnay also be seen from Figure 4.6a and 
4. 6b. 

The above values of optimum feedback gains are 
determined for a rigid airplane. It may be of interest to 
study the response of a flexible airplane with the realistic 
actuator with optimum elevon deflection feedback gain as 
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determined above. Since the model of the flexible airplane is 
not available for carrying out such a study* it was decided 
to reduce arbitrarily the efficiency of the el e von by 20% 
Cas though effected by the chordwise flexure of the control 
surface:) without changing the state matrices. Figure 4.7a and 
4.7b show the airplane response with a reduced efficiency 
eleven. The normal acceleration response of the **flexible*' 
airplane shows slight waviriess* though the settling time 
compared to the response of the airplane with ideal actuator 
remains the same Ci.e. 3 seconds^. However* the pitch rate 
response CFigure 4.7b:> shows reduction in damping and ibe 
response is osci ilatory. This is evident for the optimal 
eleven position feedback gain was determined on the basis of 
normal acceleration response. 

A study is made to observe the effect of the actuator 
velocity feedback on the airplane response. As is seen from 
Figure 3. 33b the actuator response showed marked improvement 
when the velocity feedback was introduced. However* no 
significant improvement is observed in the airplane response 
as may be seen from Figure 4. 8a and 4- 8b. The overshoot 
of the normal acceleration as well as the pitch rate of the 
ai r pi ane wl th an addi ti onal actuator vel oci ty f eedback 1 oop 
increases compared to that without the velocity feedback. 
The settling time is however the same. 
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CHAPTER 5 CONCLUSIONS 


Based on the study carried out on the electro hydraulic 
eleven actuation system for the fighter aircraft. the 
following conclusions jnay be arrived at: 

The direct drive valve may be represented by a first 
order lag transfer function for linear analysis purposes. 

Increasing the valve flow gain about the nominal value 
tends to destabilise the E H actuator. The damped frequency 
remains more or less the same. The bandwidth increases. 

Changes in valve flow pressure coefficient do not show 
marked effect on any dynamical characteristics of the E H 
actuator. There is a very small increase in damping ratio as 
the value increases. 

E H actuator displays instability at low values C< GSfi 
of nominal valued? of effective bulk modulus of the hydraulic 
fluid. Increasing the value of bulk modulus Increases the 
natural frequency and the damping ratio of the system. 

Coulomb friction between the piston and the cylinder 
surface displays a distinct damping effect on the otherwise 
lightly damped system, totally suppressing the oscillations 
at the nominal value of the level of the friction. The effect 
of the Coulomb friction is more pronounced for smaller 
oscillations of the actuator piston. For smaller command 
input values, equivalent viscous friction may be used to 
represent the Coulomb friction to facilitate linear analysis. 
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A rigidly mounted E H actuator with no load has a first 
order lag response. The bandwidth is 80 rad/s. However* for 
finite nominal value of the mounting stiffness, the response 
has an overshoot and the rise time is shorter. 

Effect of variation of horn stiffness about the I'lominal 
value is visible only for very low values C10% of nominal 
valued. Damping reduces for low values of horn stiffness. 

The inertia of the control surface as a load on the E H 
actuator system brings about instability. Introducing dynaraic 
pressure feedback through a washout circuit stabilises the 
system, and yet the best step response of the system 
obtainable in this configuration is very lightly damped 
C damping ratio of O. 113;> which is far from the 
first-order -lag response expected from a well designed 
actuator . 

Introduction of actuator velocity feedback improves the 
poor dynamical characteristics of the elevon actuator system 
markedly, suppressiiig the lightly damped oscillations. 
However, its contribution to improvement of the overall 
dynamical characteristics of the airplane, as the elevon 
actuation system is integrated into the flight control loop, 
is insignificant. 

In order that the sensitivity of the deflection of the 
control surface to applied loads Cexternally applied or 
internally generated!) is minimised, the impedance of the 
elevon actuation system should be as large as possible. 

The transfer function Cbetween command input voltaae and 
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Lhe eievon deliecLionD of Lhe elevon act-uat-lon svst-a-ni, and 
its compliance are related through a transfer function 
between the force developed by the actuator for a command 
input voltage when the piston rod of the actuator is locked. 

The normal acceleration response of the airplane with 
the realistic actuator matches quite well with that of the 
airplane with an ideal actuator when the position feedback 
gain is optimised. The pitch rate response of the airplane is 
not satisfactory even for the ideal actuator. Actuator with 
the optinvLsed position feedback gain achieves nearly equally 
good performance although the gain value has been obtained to 
achieve a good nornval acceleration feedback jnatch. 

Replacing actuator position feedback with elevon 
position feedback in the flight control loop does not show 
any significant improvement in the dynamical characteristics 
of the aircraft. 

Reducing the effectiveness of the control surface due to 
possible aeroelastic effects results in some deterioration of 
the response characteristics of the airplane, but apparently 
not significantly. 

Within the limitations of the present study C confined to 
the short period dynamics of the aircraft considered!), the 
proposed E H elevon actuation system appears to result in 
overall dynamical characteristics of the airplane not 
significantly different from those when an ideal actuator 
Cpure qain!) is installed. 

uEnTK^L L>-RAm 

\ » T., <A SiPUR 

A IIWW 
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recommendati ON5;: for further worki 

The present work may be extended further as follows: 

The inertia ot the hydraulic fluid introduces a time 
delay in the system. The E H actuator n\ay be investigated 
further taking into account such time delay. 

Effect of Coulomb friction and saturation limits 
C brought about by mechanical stops to the valve and actuator 
piston displ acejnentl) on the eleven actuation system may be 
studied by nonlinear simulation. 

The effect of eleven dampina coefficients C, and on 

a o 

the dynamical char acter i sti cs of the eleven actuation system 
have not been included in the present study for want of data. 
Their effects rnay be Investigated. 

The dynamic pressure feedback circuit parameters and the 
tuning gai n have been determined by trial and error in the 
present work. A systematic optimisation ri\ay be carried out to 
determine these values. A controller different from the one 
chosen for dynamic pressure feedback loop n\ay be considered 
for a more refined study of the system. Sindlarly the tuning 
gain may be replaced by a controller other than a 
proportional controller used in the present study. 

To study the full impact of i ncorporati ng a realistic 
actuator in the flight control loop, the total airplane model 
C including both the longitudinal and the lateral dyi'iamicsl? 
may be s 1 mui a ted . An ‘’el as t i c ” air pi ane may al so be 


consider ed. 
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The present study confines itself to unit step input to 
study the short period response of the airplane. Airplane 
response to various other inputs, such as ramp inputs or 
spectral inputs, may also be investigated. 
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APPENDIX A 


A1 - 


5 or der 


DDV Transfer Function 


37.87s 4. 1.3^10 s 4 1.08^0 


8 


8. 5a»10®s^ •+ 1.785*10 + 6. SSs"^ + 1.94*10%^ + 2. 06*10' s -»• 4.25*10 


.-3_4 


4 2 


ROOTS of E)eriou\inator : 


-a. oss^i o: 


-1.008^10 a. sss^io i 

~1 . - 4838^1 
- 8 . 3053 ^ 10 '^ 


A8. Transfer Function of Electro Hydraulic Actuator Cat no loadj 

no equi val ent vi scoci tyZ) 


—13 3 “10 8 —7 9 

1.4811^10 s 4 3.4985^10 s 4 1,3411^10 s 4 1.3805^10 


+ 2.55*10^3^ + 4.3865*10^3^ + 6.125*10® 


EIGENVALUES : -5.0068 ± 649.34 i 

-144. 61 
-lOO. 45 
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A3. Transfer Function of Electro Hydraulic Actuator Cno load; 

wi th equi val ent vi scoci ty) 

5. 6843^1 + 1.4901»elO~® s + 1.3862»elO^ 


s"^ + 4. 5846»el0‘^s’^ + 4. 9053*1 O^s"^ + 1. 043*1 0®s + 6.125*10® 


EIGENVALUES: -109.87 ± 641.52 i 

-119.36 ±14.58 1 


A4. Transfer Function of Electro Hydraulic actuator with flexibl' 
mounting Cno load ; with equivalent viscosity^ 


1.1369*10 + 1.455*10 + 1.341*10® 


+ 4. 43*10®s® + 9. 1166*10'^S® + 1. 0087*1 O^s + 5.9286*10® 


EIGENVALUES: 


-66.58 ± 122.66 1 
-154.98 ± 80.09 i 


A5 Transfer Function of the Actuator Eleven Servo System 

C without dynamic pressure feedback!) 


6 

s 


+ 443.13 


5 

s 


3 . 395 s 


+ 7. 9099*1 O^s'^ 


4.3675 s -«■ 1.984*10^^ 

+ 1. 9742*1 0®s® + 1.571 3*1 0®s® + 
1.571 3*1 0®s + 4.566*10^® 


-100.72 ± 854.44 1 
—234 . 21 
-28. 35 

+10.44 ± 28.63 i 


EIGENVALUES; 
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A6 Transfer Function of the Actuator Eleven Servo System 
Cwith dy!"jainic pressure feedback!) 


4.0x10 •+ 4. 


+ 4.43xlO‘^s® 


-6 6 

09^10 s + 8. 


R S 

+ 8. IG^lO^s '♦• 


3038x10 + 2. 84s^ + 3. Els"^ + l.QBxlO^^s 

-t- 6.6134X10^^ 

1.85xlO®s^ + 1.93xlO^®s^ + l.lQxiO^^s^ + 

4 . 74x1 + 1 . 522X1 ^ 


EIGENVALUES: -100.77 ± 855.28 i 

-118. 99 ± 102. 72 i 

“O. 32 

-1 . 81 ± 15. 92 i 
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APPENDIX B 
SYSTEM CONSTANTS 


1 

SYMBOL 

DESCRIPTION 

UNIT 

VALUE 

A 

P 

Actuator Area 


3. 142 

ft 

Fluid Bulk Modul us 

1 

lOOOOO 

F 

cmax 

Maximum value of Coulomb friction 

lb 

1 

80 

I 

El e von I ner t i a 

i 2 

1 k a m 

[ 

3. 35 

Ka 

Ac t ua t or Loop Ampl i f i er Gai n 

1 

vdc / vdc 

1 12. 5 

1 

kb 

Moun ting s t i f f ness 

kg /cm 

4 

1 4.91x10 

k 

c 

Valve flow pressure coefficient 

1 in /s 
j psi 

0. 0006535 

kia 

demodui ator gai n 

1 vdc /vac 

3. 09 

kh 

Horn stiffness 

1 ka/cm 

[ 

14. 60x1 O"^ 

k 

a 

Hinge moment per unit Angle of 

1 Nm/deq 

83. 26 

Attack- 




Hinge moment per unit eleven 

Nnv'deg 

47. 68 

deflection 



kda 

ivdt gain 

vae/i n 

1 . 43 

k 

9 

Val ve f 1 ow gai n 

i ^ , 

1 n z' s 
in 

704 

1 

Horn Ar*m Length 

m 

0. 136 

m 

P 

Piston mass 

lb 

in/s^ 

1.18 

r V 

Valve spool radius 

in 

0. 1883 

V 

Total Cylinder Volume 

in 

16. 4 
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FIG. B. 4 FREQUENCY RESPONSE OF DDV 
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FIG.2.5 VALVE-PISTON COMBINATION 
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FIG- 2. 8A ROOT LOCUS OF E H ACTUATOR; PARAMETER i VALVE FLO’ 
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FIG. 2. 8B FREQUENCY RESPONSE OF E H ACTUATOR; PARAMETER : VALVE FLOW GAIN 
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FIG. 2.9 ROOT LOCUS OF E H ACTUATOR? PARAMETER* VALVE FLOW PRESSURE COEFFIC 
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FIG. 2.1 OB FREQUENCY RESPONSE OF E H ACTUATOR; PARAMETERS BULK MODULUS 
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TIME RESPONSE OF E H ACTUATOR: VARYING STEP VOLTAGE INPUTS 



aa 



FIG, 2.15 TIME RESPONSE OF E H ACTUATORS VARYING COULOMB FRICTION MAGNITUDE 
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FIG. 2-16A ROOT LOCUS OF E H ACTUATOR: PARAMETER: MOUNTING STIFFNESS 
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FIG- 2 - 16 D FREQUENCY RESPONSE OF E H ACTUATOR; PARAMETER: MOUNTING STIFFNESS 
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FIG.2-18 sign convention 



FIG. 2. IS ROOT LOCUS OF E H ACTUATOR; PARAMETER •* MASS OF PISTON AND ELEVON 
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FIG. 2-20' ACTUATOR -E LEVON SERVO (with dynamic pressure feedback) 
















FIG- 3-31 ROOT LOCUS OF ELEVON ACTUATIOM SYSTEM (I DETERMI NATI OH 
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FIG. 2. 22 <x STEP RESPONSE OF ELEVON ACTUATION SYSTEM CWITH DPF) 
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FIG. 2. aSb STEP RESPONSE OF ELEVON ACTUATION SYS f EM CWITH VELOCITY FEEDBACIO 




FIG. 2. 23A ROOT LOCUS OF ACTUATOR-ELEVON SERVO SYSTEM; PARAMETER i hORM STIFFNESS 
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FIG. 2. 23B FJ^EQUENCY RESPONSE OF ELEVON ACTUATION SYSTEM; PARAMETER j HORN STIFFNESS 
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FIG. 3.2 IMPEDANCE CHARACTERISTICS OF ELEVON ACTUATION SYSTEM 
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FIG. 3.3 IMPEDANCE CHARACTERISTICS OF ELEVON ACTUATION SYSTEM 




^ 8KaK^rvKgPAp (Xps + I) 

Vo (s) *0 (TS+l)(T:ps+1)[(V+4P/!?p/Kc)s+43Kcl+8 KoKi rv Kg P Kpt 

















Exl. load 



FIG.Kb.SIMPLIFIED REPRESENTATION OF TRANSFER FUNCTION / 
COMPLIANCE OF ELEVON ACTUATION SYSTEM 
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FIG. 3.5 COMPARISON OF TRANSFER FUNCTION AND COMPLIANCE OF ELEVON ACTUATION SYSTEM 
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FI <3. 4.3 ROOT LOCUS OF AIRPLANE C DETERKt NATI ON OF TUNING GAIN) 
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FIG. 4. 3A STEP RESPONSE OF AIRPLANE? PARAMETER : TUNING GAIN 
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FIG. 4. 3B STEP RESPONSE OF AIRPLANE; PARAMETER i TLfNING GAIN 
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FIG. 4. 4A STEP RESPONSE OF AIRPLANE — REALISTIC VS IDEAL ACTUATOR MODEL 





FIG. 4. 4B STEP RESPONSE OF AIRPLANE — . REALISTIC VS IDEAL ACTUATOR MODEL 




FIG. 4. 5A ROOT LOCUS OF AIRPLANE? PARAMETER: ELEVON POSITION FEEDBACK GAIN 





•hH 


C 






X 




X 





rX 











in 


s 




X 

in 




u 





if 








n 




< 




If 





CQ 





n 




P 




p«t}Cn 

C) 

r^» 




W 

W 

u. 




rd 

rj 






Q 


(d 









,h(Q 




o 

1 — i 









r*i 

H 


o 

ri 

iD 

0 


\n 




■ujlUJ* 


Ch 

.0 


- 


§ 


rrj 

a 


C|«< 


ro 


P 




C 







• •Ah J 


b 

• 




? 


b 

o 

«r«t 

c-* 

/*“• 




$ 


rd 


-4-lf 

l<zm)< 




i:ia 




-pH 


- 

CO 




l— M 

• p»'4 


0 



w 


I’d 

ID 

SI 

O 

•t;5 



(A 




P— t 





X 


• p»< 

rcJ 

* 

O 




M 



iD 




! 

TIME 

H 


i*”* 

ItjfMi* 

"P-H 

3 


0 

id 

1 

0 

0 


CM 

P 

O 

Od 

O 



Mp.«l 


** 






1 

3 

1 



CM 


1 









j 







m 


& 







•ijH 



***** **’*% 








u. 









0 

'*'*••!* *• *. 










. " * . "’• 










■•-- 






z 




““ ^ 





o 







m 


Ui 

a 







... 







■“-, 


CSD 


a. 









U2 

IH 

1 

1 


1 

1 _ 

A 

_ .£ - . 

- GD 


p 

m 

• 

1 

m 

vj3 



(X3 






m 

CSD 

{ 


• 









o 









w 


MO I X 

'lowa OM 



Cih 





c 




> 




l\3l 


t 

JL 


s 

XJ. 




c > MO » .iUJI ‘IMMH OW 




0 

:rj 

. . . 

a.:i -»• 

irs ;?• ■ 

. P 0 


Sr: 

O Lp 
&' '..p 
*'>««•' . •• 

' ' ' l±)' 
II ■ 
11 


u 

.4..^ 

- 1 

1 I] 


o 

i'iff 





{Vj q t 


nj 


I..Q I'd : 



i'J 

•'ci rj;a .1 


...p 

JTj 

lU '"iri 1 

‘l 

c:i 


Cl Cl 

i • 

I'd 


'E +.4 llJ • ‘ 

I 

r 




i \ : 


..p 


r 1 

rd 

w 

0 ff: , 

» 

i • 

0 

" r***! 

«i”< Ip *, 

i **, 

-O' 

|nn.>ii 

•i^iujt * i«"i 

V • 

BimiJ 

csj 

« /««J 4 *^ 

. I ’.. 

• *‘v 


. . .flJ 

. M ■T' 



SI.! 

0 C') 



0 

: '■■•1 i 


fnMni 

A 

* \ * .. 

3 


i 



1 * |«.rl 

-- 

'x ,* 

. . i . . 

3 



• **'*■. 


< s / {} 1> ,A ) a X M H 1 1 : > I f < t 


FIG. 4. 6B STEP RESPONSE OF AIRPLANE FOR OPTIMAL ACTUATOR POSITION FEEDBACK GAIN 
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FIG. 4, 7A STEP RESPONSE OF AIRPLANE FOR REDUCED EFFICIENCY OF THE ELEVON 





FIG. 4. 7B STEP RESPONSE OF AIRPLANE FOR REDUCED EFFICIENCY OF THE ELEVON 
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8a STTEP RESPONSE OF AIRPLANE CWITH ACTUATOR VELOCITY FEEDBACK! 



llll<lSd JL f 



CD 


CS3 




LD 

»j;3 


iN- 

cc 



CD 

CD 

(D 

m 


CD 



CD 

CD 



1 

CD 

I 

I 

( 


CD 

{ 

f 


® 

{ 

i 



* 

1 

f" s 



lyjf 

11:) 

'.1 

r .1 




FIG. 4* 8b STEP RESPONSE OF AIRPLANE CWITH ACTUATOR VELOCITY FEEDBACJO 




